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Abstract 

It is shown that the arguments in the reply of Z.-D. Zhang defend- 
ing his conjectures are invaUd. His conjectures have been thoroughly 
disproved. 

After all the discussion about his paper [H El El HI [5] , Zhang seems to have 
only one real issue left with the two comments [21 [5] , still wrongly believing 
[6] that the free energy of the three-dimensional Ising model is not analytic 
at /3 = 1/(/cb^) = 0, H = 0. His further arguments are irrelevant or dealt 
with adequately in [21 [5] . 
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His objection that Gallavotti and Miracle-Sole set /3 = 1 is not valid. 
One often uses dimensionless parameters Ki = pji, {i = 1,2,3), h = f3H. 
Equivalently, one can absorb the (3 into the coupling constants, setting (3=1, 
so that Ki = Ji, h = H. Infinite temperature is then the limit Ki ^ 0, h —>■ 0, 
such that all ratios are kept fixed0 

His point on [7] is also not vahd. Writing z = exp(— 2/5if), see eq. (23) 
of [7], and keeping (311 fixed in the limit /5 — > 0, the Ising model partition 
function on an arbitrary lattice with N sites becomes Z = (z^^"^ + z~^^'^)^ , so 
that all infinite-temperature zeros of Z occur at z = —1, i.e. for the purely 
imaginary magnetic field = ±i7^kBT/2 = ±ioo. There is no T = cxd 

singularity at H = 0, z = Iq 
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^The reduced free energy per site (3f is often rewritten (3f = (j){{Ki}, h) = 0({/3Ji}, (3H) 
with some ftmction <j). Setting /3 = 1 is no loss of generality, as one can easily restore the 
/3-dependence by the replacements Ji f3Ji, H — fiH , f (3f . 

^In Zhang's paper |JL,, = 0, z = 1, Z = 2^ for any lattice with N sites at T = oo, 
which is far from z = —1. Hence, the infinite-system dimensionless free energy /?/ is 
analytic at /3 = also by the general theory of Yang and Lee [5] . 
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Added Comments 



It has saddened my heart to see 'arXiv:0812.0194V3 and to see that Zhang is 
still not convinced that his conjectured results are in error. For his benefit 
and the benefit of others who might be confused by his remarks I am adding 
the following to show why his newer statements are also in error. 

The probability interpretation of statistical mechanics involves Z and 
thus Therefore, that / — > cxo for the free energy per site as T ^ cxo is 
of no physical significance. Series and analyticity determinations must be 
starting from the reduced free energy (5j = f/k^T, not /, near T = oo. 

Also, as long as (the real part of) / is negative, (which is easily checked 
for the Ising model at arbitrary temperature), Z = exp{—Nl3f) blows up and 
l/Z— >-OasA^— >-oo. So limAT^oo l/Z = occurs for all finite temperatures 
even at the critical point where the Yang-Lee zeroes pinch the real axis. 
It is a misinterpretation of the Yang-Lee papers to study zeroes of l/Z. 
The Boltzmann probability is given by exp{—f3S)/Z, with Z a polynomial 
in z = exp{—2j3H), H being the scaled magnetic field. This probability can 
only show anomalous behavior if zeroes of Z approach as ^ cxd, as is 
clearly explained by Yang and Lee in their two famous 1952 papers. 

Finally, there are many papers proving analyticity of the reduced free 
energy in 1/T around T = oo, for both classical and quantum models that 
generalize the Ising model. Such proofs have been published by many groups 
in Europe, Russia, Japan and America. I have only quoted a few mostly 
Ising references in my comment. The exactness of the series is rigorously 
proved many times. It is also backed up by numerical work, as Fade analysis 
of the series expansion is in excellent agreement with other numerical work, 
such as Monte Carlo calculations. There can be no such evidence backing up 
Zhang's work, as it violates rigorously established theorems. 

Added Comments 2 

As Zhang is stacking up further errors in his fourth version a rXi v : 812.01 941 ^4 
and repeats his errors to the followers of his blog in China, I feel that I have 
little choice than adding another reply. 

In the first place, under "solving the three-dimensional Ising model" is 
understood the calculation of some basic thermodynamic quantities using 
standard equilibrium statistical mechanics starting from H in [1], i.e. the 
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usual Ising interaction energy. Any talk about time dependence, ergodicity, 
relativity, quantum mechanics, black holes, etc., is a distraction from the 
wcll-dcfincd problem and changing it to a different problem. 

Standard equilibrium statistical mechanics states that the free energy per 
site / is calculated using 

M=j^^'' = -j^iogZ, (1) 

with Z the partition function and N the number of sites. It is /9/ which is 
calculated first and it is /?/ for which series expansions are done. That / 
becomes oo at T = oo is of no significance, as is obvious from one of the 
basic formulae in thermodynamics, F = U — TS. This says that at T = 
we have F = U the (internal) energy, while at T = oo we find F/T = —S, 
i.e. minus the entropy. This explains again that F/T, or equivalently (3F, is 
preferred over F at high temperatures. 

Next, consider eq. (1) in Zhang's paper [1], and define 

C = max{\J\,\J'\,\J"\}, (2) 

the absolute maximum of the coupling constants. Then it is easily checked 
that 

1^1 < 3NC, \Z\ < 2^e^^^, log |Z| < log 2 + 3C, (3) 

as the number of terms in Z is 2^ . The only way to get singularities in (5f 
is the occurrence of zeroes of Z, as explained excellently by Yang and Lee, 
whom Zhang chooses to misquote. Zhang's work is in stark contradiction 
with the work of Yang and Lee result. 

Finally, using the logic of his conjecture we can solve almost any problem 
in statistical mechanics, but with incorrect results. The idea of going to 
higher dimension is a generalization of the fundamental theorem of calculus 
JI^ f'{x)dx = f{b) — f{a), with examples the theorems of Green, Stokes, 
Gauss, and the Wess-Zumino term. For the "guess" to work, Zhang must 
show that the integrand in four dimensions is a derivative of the integrand in 
3D (or a discrete version with integrand replaced by summand and derivative 
by some kind of difference). In [1] I cannot find that Zhang understands this 
and, as the conjecture is rigorously proved to be false, it makes no sense to 
further discuss this matter here. 
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Added Comments 3 



Most proofs of the analyticity of free energies and correlation functions use 
linear correlation identities of Scliwinger-Dyson type, known under such 
names as the BBGKY hierarchy, Mayer-Montroll or Kirkwood-Salzburg 
equations. If Zhang had read [9j in detail and understood, he would have 
found more than one way to arrive at a proof of analyticity of f3f at = 0. 
He would also have found [lOj (ref. [12] cited in 0). Let me try to explain 
the proof in more down-to-earth terms using an identity of Suzuki [111 [12] , 
restricted to the isotropic Ising model on a simple cubic lattice with periodic 
boundary conditions and of arbitrary size, i.e. 

n^..) = ^E(fn-.0tanhf/3J E ^i))^ (4) 



i=l k=l / nn of jfc 

where ji , . . . , j„ are the labels of n spins and I runs through the labels of 
the six spins that are nearest neighbors of aj^,. The isotropy assumption, i.e. 
Ji = J2 = J3 = J, is made to simplify the argument, but can be easily lifted. 
Averaging over k has been added in (jl]) , so that all spins are treated equally. 
Next we use 

tanh f/JJ^o-H = aij^f^i + (^sj^^h^^h^^h + (^^J^^h^^h^^h^^k^^h^ (5) 

^ 1=1 ^ (6) (20) (6) 

where the sums are over the 6, 20, or 6 choices of choosing 1, 3, or 5 spins 
from the given cxi, . . . , ae. It is easy to check that the coefficients are 

t(l + 16^2 + 46t^ + 16^6 + ^8) _ _2t3 



16t^ 

The poles of the a, are at t = ±i, t = ±(\/2 ± l)i, and t = ±(v^ ± 2)i. It 
can also be verified, e.g. expanding the in partial fractions, that the series 
expansions of the Oj in terms of the odd powers of t alternate in sign and 
converge absolutely as long as \PJ\ < arctan(2 — v^) = 7r/12. 

The system of equations (jl])-© can be viewed as a linear operator on 
the vector space of all correlation functions of the 3d Ising model. It is easy 
to estimate the norm of this operator, from the 32?t, terms in the right-hand 
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side (RHS) of (jl]) after applying using the alternating sign property of 
the Oj's. It follows that we only need to study 



oai + 20a3 + 605 = — 7- (7) 

for purely imaginary t to find the desired upper bound r for the norm. We 
then have that the RHS of (jl]) is bounded by rM, where M = max \ {a ■ ■ ■ a)\ 
with the maximum taken over all 32n pair correlations in the RHS. Obviously, 
M < 1 if /5 > and real. We must stress that the bound is also valid if all 
the Oj's are replaced by the power series in t obtained from the power series 
of the Oj's replacing each term by its absolute value. 

According to the above we set r equal to the absolute value of ([7]) for 
imaginary t. We can then show that r < 1 for 

\t\ <{V3-V2){V2 - 1)= 0.131652497 ■ ■ ■ , or 

\PJ\ < arctan[(V3 - V2)(V2 - 1)] = 0.130899693 ■ ■ • . (8) 

To prove analyticity of /?/ in terms of /? at /3 = it suffices to study the 
internal energy per site or the nearest-neighbor pair correlation function, as 

u = -^^p- = -3 J (croooO"ioo)- (9) 



We apply (jl]) to (croooO"ioo); then we apply (jl]) on each of the new correlations, 
and we keep repeating this process. We may now and then encounter a 
correlation with zero a factors (1) = 1, at which the process ends. As each 
other correlation vanishes with a power of t, we then generate the high- 
temperature power series to higher and higher orders, for arbitrary given 
size of the system. The series is absolutely convergent as the sum of the 
absolute values is bounded by < when (jH]) holds. We can at first 
assume that /3 > and real. But from the absolute convergence we also 
obtain a finite radius of convergence in the complex t and f3 planes. 

Increasing the size of the system, we conclude that more and more terms 
become independent of the size, whereas the remainder rapidly tends to zero. 
Therefore, the series converges to the thermodynamic limit, and we have once 
again proved that the free energy and all correlation functions are analytic 
in t and in (3, as long as (jH]) holds. This bound is a rough estimate and 
much better ones have been given in the literature. Furthermore, adding a 
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small magnetic field H and generalizing the steps in the above, we can also 
conclude that all correlation functions are finite for small enough f3 and H, 
so that there are no Yang-Lee zeroes near the H = axis for small j3. 

In conclusion once more, papers [U [3], [6] contain serious errors and are 
beyond repair. 
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